IFUP-TH 59/92 



Finite size scaling in CP models 



Paolo Rossi and Ettore Vicari 
Dipartimento di Fisica delVUniversitd and Istituto Nazionale di Fisica Nucleare, 1-56126 Pisa, 
Italy 



Finite size effects in Euclidean CP^ -1 models with periodic boundary conditions 
are investigated by means of the 1/N expansion and by Monte Carlo simulations. 
Analytic and numerical results for magnetic susceptibility and correlation length are 
compared and found to agree for small volumes. For large volumes a discrepancy is 
found and explained as an effect of the physical bound state extension. 

The leading order finite size effects on the Abelian string tension are computed 
and compared with simulations finding agreement. Finite size dependence of topo- 
logical quantities is also discussed. 



PACS numbers: 11.15 Ha, 11.15 Pg, 75.10 Hk 



Typeset Using REVTEX 



1 



I. INTRODUCTION 



The study of the finite size effects has become an important tool in the theoretical and 
numerical analysis of critical systems and lattice field theories. 

Finite size scaling may be used in order to improve the extrapolation of finite lattice 
results to the infinite volume limit relevant to field-theoretical predictions. More naively, 
knowledge of the volume dependence of lattice results in sufficiently large volumes allows to 
choose lattices such that finite size effects may be consistently smaller than the intrinsical 
error of numerical simulations. 

Within our program of studying analytically and numerically two-dimensional CP^ -1 
models as a prototype for lattice QCD, we decided to explore their finite-volume behavior 
with special emphasis on two peculiar properties that our model share with QCD: existence 
of a linear confining potential and non-trivial topological properties. 

CP^" 1 models are 1/N expandable: we can therefore, for sufficiently large N, make 
theoretical predictions on the size and pattern of finite lattice effects. This chance, in 
conjunction with the opportunity of performing high-statistics Monte Carlo simulations for 
intermediate values of iV with an acceptable computer effort, offers the possibility of a 
systematic study of finite size effects. 

This papers is organized as follows: 

In Sec. [II]we review the theory of finite size scaling in the context of the 1/N expansion. 

In Sec. |TTT| we discuss the lattice formulation adopted, introduce the physical observ- 
ables (magnetic susceptibility and correlation length) and evaluate the 0(1) and 0(1/ N) 
contributions to their finite size functions. 

In Sec. |T^ we perform the same analysis for the leading order contribution to the Abelian 
string tension. 

In Sec. |V| we present the results of numerical Monte Carlo simulations for selected values 
of N, from N = 2 to iV = 100, and compare them to our theoretical predictions, explaining 
the discrepancies, when present. 

In Sec. |VI] we draw some (hopefully general) conclusions. 

II. FINITE SIZE SCALING AND 1/N EXPANSION 

Finite size scaling in the 1/N expansion has been introduced and discussed in Ref. [[!]. 
We recall here the main results adapting them to the specific context of CP^ -1 models. Any 
coordinate- independent physical quantity will in general depend on four different parameters: 

Q = Q(g,a,L,N) (1) 

where g is the coupling, L d is the physical volume in d dimensions and a the lattice spacing. 
However in the critical region g — > g c and in the infinite volume limit all dependence on a 
and g may be absorbed in the dependence on the correlation length 
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The finite size scaling relation amounts to the statement that, when £ — > oo while keeping 
L/£ finite, we should get 

Q(g,a,L,N) f{Q) 



Q(g,a,oo,N) 
The 1/N expandability implies that, assuming 



fW(L/£,N) . (3) 



Q(g, a, L, N) = Q (Ng, a, L) + ^(Ng, a ,L)+o(-^) ( ! ) 



Z(g 7 N) = £ (Ng) + h 1 (Ng) + o(±-\ , (5) 



and 

we may expand the finite size functions in the form: 

/W)(L/£,AT) = / ^(L/0 + ^A (Q) (VO + 0(^) , (6) 
and obtain the representations: 

®f N g ' a ' L \ - (7) 

Qo(Ng,a,oo) 

Qi(Ng, a, L) Q^Ng, a, oo) | L /q (Q) (L/£ ) 6 _^ A (Q) (^o) (g) 



III. LATTICE FORMULATION 

Our choice of a lattice action for CP^ 1 models is @: 

^ ] (^n+fiZ n \ n ^ ~\~ z n z n j r ^ L \ n ^ 2^ , g , (9) 

where z n is an iV-component complex scalar field, constrained by the condition z n z n = 1 
and X nifl is a U(l) gauge field satisfying A n>Ai A njAt = 1 . 

Observable physical quantities must be gauge invariant; we therefore focus on the corre- 
lations of the composite operator 

Pij(x) = Zi(x)Zj(x) , (10) 

and specifically on the two-point function 

G P (x) = (TtP(x)P(0)) c . (11) 

Starting from the continuum definition of magnetic susceptibility 
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Xm = J d 2 x G P {x) 
and (1/N expandable) correlation length M 



e 



Jd 2 x \x 2 G P {x) 
J d 2 x Gp(x) 



we can introduce their counterparts on finite lattices with periodic boundary condit 
us first define the finite lattice Fourier transform of the correlation function 



G P {k-L) = -^(TrP(m)P(n)) c exp 



2m 



(m — n) ■ k 



defined on integer values of the components of k. We then define 

Xm(L) = G P (0, 0; L) — ► Xm , 

L— >oo 



em 



i 



4 sin 2 tt/L 



Gj(0,0;L) 
G P (1,0;L) 



L— >oo 



e ■ 



Within the 1/N expansion the following relationship holds 



G P (k;L) = r 2 

where, in the infinite lattice notations: 

d 2 q 



^)(k;L) + ^\k;L) + 0(^ 



\N 2 



A ^ = / (ar) > 



[q 2 + ml] 



(k + q) + ml 



A^(k) 



(Jlq A {a) (q)V {a) (q,k)+ I ^- 2 A w ( q )V {x) (q,k) 



(2n) 2 

d . _i s /" d 2 q 



(27T) 



+ \^B^ A ^i '(ar). 

The auxiliary quantities we have introduced are the effective vertices 

* _i / , \ /• d 2 g 2E^cosg^ r rf 2 g 4 E M sin 2 + fc M /2) 
(A)l ^ " ■/ (2-) 2 ? + m 2 J (2n) 2 [ga + ^ [ r „~^ + „ 



+ fc) +m„ 



V(a)(^,P) 



d 2 q 



(2ny 



[q + ml] 2 



(q+p) + m( 



d 2 q 



(q + k) + ml 
1 



+ p 



~P 



^ 11 1 [q + ml] (q + p) + ml (q + k) + ml 



(q + p + k) + m 2 , 



V(x)(k, P ) 



d 2 q 



4E M sin 2 (g / , + ^/2) 



+ 



(2?r)2 [q + mlf 
d 2 q 



+ p) + m 2 



(q + k) + ml 



+ p 



~P 



4 £ sin(g M + p M /2) sin(g A1 + fc M + 



( 27r ) 2 [g + mo] (g + p) + (g + k) + m 2 , 



(g + p + A;) +mg 



d 2 q 



2 E M cos g^ 



(27T) 



[g + mg 



2l2 



(22) 



The conversion to finite lattice is obtained by the replacements 



d 2 q 
(2vr) 2 



L 2 



E , 



(23) 



and on periodic (toroidal) lattices the zero modes of the inverse propagator A^(k), if 
nonvanishing, must be explictly removed. Moreover 



777 



L i 



where in order to keep Ng fixed we must relate mi to mo by the lattice gap equation 



-E 

T2 

9^ = 1 



L 2 g 2 + m 2 



1 



d 2 g 



A" 



(2tt) 2 g 2 + m 2 2nl+m 2 /A \ 1 + 



(24) 



(25) 



This equation has been analyzed in detail in Ref. M for m — * 0, i.e. in the scaling region, 
where 



d 2 q 



1 



1 , 32 
in ■ 



(27r) 2 q 2 + m\ m -.o 47r mj 
In the scaling region one could establish the relationship 

z = m -L = z c exp 



(26) 



(27) 



where zl = ttilL, z c = 4.163948..., and in the region zl < 27r the function u> may be defined 
by 



d„ 



Air 



+ 47r£(-l) n *f n d n+ i , 



n=l 
oo 



(2vT 



E 



n > 1 



(28) 



In the definition of d n the sum must be performed excluding the term (m 1; m 2 ) = (0, 0). The 
function Zq(zl) is monotonic and invertible. We could therefore make use of the auxiliary 
variable zl in the computations and replace it by zq only at the end. 
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In practice, in order to take into account at least some of the scaling violation effects due 
to irrelevant operators, we found it always convenient to solve explicitly and numerically 
Eq. ( p5|) on each finite L x L lattice for assigned value of Ng. 

Let us derive the large- N results: 



(3 2 Xm , (L)=A7'(0,0;L) = ^£ ' 



(3 2 Xm,0 



d 2 q 



L 2 ^ (q 2 + m 2 L ) 2 
1 1 



d{3 



(27r) 2 (g2 + m 2)2 mo ^ 4?rm 2 



,2 ' 



(29) 



and 



4 sin 2 tt/L 



£ 2 - 



lim 



L— >oo 



Ag)(0,0;L) 
^(1, 0;L) 
1 



- 1 



m -*0 Qrrin 



(30) 



The corresponding large- iV finite size functions are obtained by the relations 



£o(&oo) 



/o (0 (£/£o) 



(31) 



for "~ > °° while keeping L/^ constant. 

The infinite volume limit quantities can be obtained numerically by performing the cal- 
culation on very large lattices and checking the convergence of the results. For large enough 
lattices an extrapolation based on the formula 



Q(L) — -> Q(oo) + 



L 2 



for L > 1 



(32) 



gives safe results even if it does not keep into account the existing logarithmic corrections to 
the L~ 2 term. Furthermore, we must in principle worry about the effects due to irrelevant 
operators that may spoil the universal behavior by introducing a separate dependence on 
£o~ 2 . Working at very large values of £o takes partially into account this problem. Further 
improvement may be obtained by considering different values of £o and linearly extrapolating 
in £^ 2 . For £ large enough, this extrapolation, albeit purely phenomenological because all 
existing logarithmic dependence on £ is not included, is powerful enough to predict results 
that are stable under iteration of the procedure (within the numerical precision required). 

An accurate evaluation of the functions /q and /q is presented in Figs. [I] and § 
respectively. The behaviors are quite similar and are characterized by the existence of two 
different regimes, corresponding to different physical situations. As long as L/£q is small 
{L < 5£o say ), we are exploring the small volume behavior of the system. Because of 
asymptotic freedom, this is an essentially perturbative, small coupling regime, and actually 
the finite size functions might be quite accurately evaluated within standard finite volume 
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perturbation theory [|]-||, with the only limitation of ignoring the absolute scale of the 
abscissa, since perturbation theory does not allow for the computation of the ratio between 
the mass gap and the A parameter which relates the quantity AL appearing in perturbation 
theory to the finite size variable L/£ . When L/£ becomes large (L > 10£ say ) we start 
to explore large-distance effects, that are dominated by the physical mass poles. Due to 
the two-particle nature of the operator (|10D , the functional dependence is characterized by 
a decrease with growing L/£o (two particles in a box effect) in contrast with the asymptotic 
growth characteristic of the corresponding single-particle operators (free-particle in a box 
effect). 

In the scaling region we can write in terms of oj(zl); indeed we have 



JO i L t,0 > n ^— 2 - ~ Z o 2 ' i 66 



which is useful for zl < 2tt. A similar formula can be obtained for fff\ 

The determination of the 1/N corrections to the finite-size functions proceeds according 
to the general formalism [|l| described in the previous section. Let us consider as an example 
the function /i ; the determination of f[ involves exactly the same steps and the details 
will not be reported here. We may deduce from Eqs. (P|, |r5] , |TTD the relationship 

f ( o X \L/^) ^(0; L) A^(0,oo) £o / «(L/£ ) £o ' 
where one can prove that 

A "<°) = 4/|W> (0)A( ° ,(,) m ( a '"» (o) ^ a -» (9) ) 

+A^(0)A w (,)^j (a^O^A^,)) , (35) 
and its finite lattice generalization holds. The function 

" h, (£/?o) = ItStM = -5? ln/ ° rf (36) 

could be extracted from the definition of /q 0^/£o)- in particular one can show that in the 
scaling region 



t _ L 2 E(q 2 + m 2 L )- 3 1 d 2 u/(dz 2x2 

Finally we can evaluate 



# (x) (V£o) — 1 - - ^ - T 2 = 1 - - w/ V "J 2 ■ (37) 

™0-0 27T[^(g2 + m 2 ) _2]2 {du/dz 2 ) 2 ^ ' 



7^ = ^^ AL) 



A^(l,0;L)4sin 2 7r/L 



Aj" 1 (0,0;L) A^L) 
A^(0,0;L) A^(1,0;L) 



(38) 
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We are now in principle ready to evaluate the 1/N contribution to the finite size functions 
of the magnetic susceptibility and of the correlation length. 

There is a further subtlety one must take into account when evaluating quantities in a 
finite lattice. The definition of H^ x ' implies that this quantity must vanish in the L/£ — > oo 
limit. Contributions from irrelevant operators may therefore become comparable to the value 
of the function itself in the large L regime. We must make sure that these contributions do 
not spoil the final result. A very simple trick amounts to replacing the definition extracted 
from Eq. (K7l) with the following 



L I r_5_£(g2 _l™2^- 2 1 2 II \f(a 2 + 



rn 



m 2 L/ 



{[Kg 2 



+ml)- 2 ] 



(39) 



Since in the scaling region 



J(q 2 + m 2 ) 

U(q 2 + ml) 



mo— >0 



2tt 



(40) 



the two definitions differ only by irrelevant contributions. However the second definition has 
the property that it goes to zero when L — > oo even outside the scaling region, and therefore 
it will not spoil the large L behavior of the function. 

The same analysis goes through in the construction of the finite size scaling function f± 
and we will not repeat it here. Let us only quote the final expression we adopted in the 
numerical calculations: 



z 2 d\nfP/dz 2 



1 



+ 



1 



1 



2 4sin 2 7r/L 4vr£ 2 (L)' 



A (Q) (1,0;L) x 



A(«)(l,0;L)-5: 



[q 2 + m 2 



(Q + l)' 



+ m\ 



A (Q) (0,0;L) — J2 



L 2 [q 2 + m\ 



, (41) 



dO i AO 



where I = (1, 0). 

We are now ready for a numerical evaluation of fi^/fo^ and f\ i> /fo' > - The results are 
presented in Figs. |3] and [| respectively. Values obtained at finite £ 0) an d therefore affected by 
the presence of non-scaling contributions, are subsequently extrapolated at infinite £ > after 
having verified that scaling volations are well parametrized by a term proportional to £^~ 2 . 
We estimate the error in our determinations to be uniform and to be everywhere contained 
within the graphical width of the line representing our results. 

We notice that the behaviors are similar but they are not quite the same. In particular 
fi^ has two zeros: the smaller one lies in the perturbative region and should manifest itself 
in the existence of a region where the finite size function at finite iV is, for iV large enough, 
independent of N. This result will be confirmed by Monte Carlo simulations in Sec. [V]. 



IV. STRING TENSION 



Another very important physical quantity in the context of CP^ 1 models (and of lattice 
gauge theories) is the string tension, that may be extracted from the asymptotic area law 
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holding for Wilson loops. Within the 1/N expansion, the string tension is depressed by a 
power of iV _1 , and therefore it is sufficient for our purpose to focus on the first non-trivial 
contribution to the rectangular Wilson loop: 

1 r r\ 2 h sin 2 ^-^ sin 2 ^ 
*W{R,T) = -±- / Sm 2 2 fe 2 A (A) (fr) . (42) 

2iV J {2ir) z sin y sin y 

Actually in order to extract the string tension it is convenient to study the Creutz ratios: 

W(R,T-1)W(R-1,T) 



X c(R,T)= In 



W(R,T)W(R-1,T-1) 
1 ^ M |(2H-l )sm |(2T 1) p 
i (27T) 2 sin^ sin? v ; 



2iV7 (2tt) 2 sinf sin f 

and the Polyakov ratios 

W(R-1,L) 



X p(R)= In 



W(i?,L) 

"2 



1 /■ difei sin %(2R- 1)— 2 A „ , 
= 2n!^ 'L$ *■*"><*■■ 0) • (44) 

In the infinite volume limit and in the scaling region, for R,T ^> £ one obtains 

Xc(i?,T)^— = ^ = a 

Xp(R)-^v ■ (45) 

On a finite lattice xc* an d Xp are obtained by replacing in Eqs. fl43]j44l) A(a)(/c) with the 
finite size propagator A^(k;L). We may focus on the square Wilson loops (R = T) and 
compute the dimensionless quantity C(R) = Xc(R, R)£, 2 - However, when studying finite 
size effects, we must take into account the existence of two different physical length scales: 
£ and R. Therefore, even in the scaling region we end up in a function of two variables; 
indeed for £ — > oo keeping fixed the ratios L/£ and R/£ we must have 

C(N,g,R,L,a) — f^(N,L/^,R/L) — > ±-fP{L/£ ,R/L) . (46) 

The same relation holds for the dimensionless quantity P(R) = Xp{R)^ 2 '- 

P(N,g,R,L,a) f^(N,L/^,R/L) ^-/ (P) (^o,i?/^) • (47) 

£^oo iV— >oo JV 

We evaluated the function /q numerically at the fixed value (3 = 0.86 such that £o ~ 16, 
where we checked that scaling was satisfied within 0.2%. The results are plotted in Fig. |5| 
for several values of L as a function of R/L. The resulting curves correspond to the sections 
at constant L/£o of the two- variable function /q^. This choice of variables is most adequate 
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if we want to focus on large Wilson loops (i.e. such that R £). Otherwise a more proper 
choice would have been the couple L/£o> R/£o- 

Notice that even in the limit L/£q, R/£,q — > oo as long as the ratio R/L is nonvanishing, 
we do get a nontrivial function of this ratio. Since at large N, in the continuum and for 
infinite volume, we have that for fc£o — > 



A(A)(*) 



2vr 



(48) 



we may exploit the condition L/£ ^ oo with R/L finite to replace the dressed propagator 
A( A )(/c; L) with a lattice representation of its pole part |7|]: 



A em (fc) — 

It is then easy to show that 

limx/^oo fP(L/£ ,R/L) 



2n 



(49) 



7T 



1 



2i2- 1 
Z 



(50) 



while the corresponding computation for the Polyakov loops leads to 

■2R - V 



lini/. c : .x pl 



ff>{L/t Q ,R/L) = 7T 



L 



(51) 



Actually one may easily prove that the vanishing of the Creutz and Polyakov ratios at 
R = (L + l)/2 holds for any finite lattice and for any acceptable choice of the propagator 
A(a), and depends only on the Abelian structure of the loop variables. 
We may study the lattice symmetries and further notice that 



fP(R) = -fP(L + l-R) 
while, defining the combination 

Xsym(R,T) = X c(R,T) - X p{R) - Xp(T) 



we find 



Hence 



Xsym(R,T) = Xsym{L + 1 - R, L + 1 - T) . 



fo ym (R) = fr(L + l-R) r7r -7T 



1 - 



2R 



L 



(52) 



(53) 



(54) 



(55) 



It is tempting to redefine the finite size scaling functions extracting the purely kinemati- 
cal factors derived above. We shall consider for definiteness the quantity (even under the 
exchange R — > L + 1 — R) 
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f P(L/Co,R/L) 



When L > ( we can study the regime L > 2R 3> £ 0) expecting our results to be largely 
independent of the ratio R/^q. Naively this might appear to imply that gQ P \L/l; , R/L) 
should not depend on R/L. However this condition is realized in the following subtler way: 

gP(L/£ , R/L) — vr + <f>(L/^(R/L) (57) 

for L/£o ^ 1 an d R/£o > 1; an d the function 0(L/£ O ) — > for L/£o — > oo. 
(p) 

The function has been numerically obtained directly from the abovementioned de- 
termination of and is plotted in Fig. ||. We extracted our best determinations of the 
functions <p an d t/> from an analysis of the same numerical data, when they were relevant 
(i.e. for R/C,o > !)■ Our results are presented in Figs. |7| and |8|, where we have chosen the 
normalization ip(l/2) = 1. Studying the Creutz ratios we found a similar behavior. 

The origin of the nontrivial albeit kinematical (i.e. not related to the propagation of the 
physical degrees of freedom) dependence in the regime R/£o > 1 can be traced qualitatively 
to the observation that, for all finite lattices, the dynamically generated massless pole in the 
propagator A^(k) does in fact disappear, and the linear confining potential is replaced, in 
the small k regime, by an effective Yukawa potential parametrized by a "mass" fi = /z(L/£o) 
becoming smaller and smaller with increasing L, and therefore approximating the linear 
behavior with better and better accuracy according to the relationship 

1 — e~ MiJ 

lim^o = R ■ (58) 

The parameter /j(L/£o) sets the scale of an effective long range lattice interaction; when 
L/£ — > oo /z(L/£o) — ► 0, and the infinite volume physics is well approximated when 
1/fi ^> L; but when L/£o is not too large the relationship L > 1/fi ^> £ may hold. 
As a consequence of this phenomenon, the measured finite lattice string tension, whose 
operational definition might be 

g { P \L/Z Q , R/L = 1/2) — vr + <f>(L/£ ) , (59) 

may turn out to be significantly different from ir even for large values of L/£q. 

In order to have a more quantitative perception of this phenomenon, we performed a 
study of the finite size dependence of the function AT^ (k, 0; L), plotted in Fig. |^ as a function 
of the dimensionless variable kC,o for several values of L/£o- Notice that the variable k£o can 
only take discrete values, and for small L/£o these values are quite distant from each other. 
In Fig. [L0] we plot the values of A^(0; L) = fi 2 ^l/2^ against the function (£ /L) 2 taken as 



an estimator of the contribution from the first nontrivial value of k. Fig. |10] shows that the 
condition A7^(0;L) <C ((, /L) 2 (corresponding to L <C 1/fi) holds for L/£ > 25, while for 
smaller L/£o significant deviations from the infinite volume behavior are to be expected, as 
already phenomenologically observed. 

We evaluated analitically the function gy'\L/^o, R/ L) in the scaling limit, assuming an 
effective Yukawa propagator, and found 



11 



where 



(/iL) 2 = 27rA (A 1 ) (0;L/6 




(61) 



The corresponding functions 



<MV£o) = t 



T 



sh(f; 



(62) 



and i/jy(R/ ' L) are plotted in Figs. [TT] and [12|, using the same conventions and notations of 
Figs. |7|, H The qualitative and quantitative agreement is extremely satisfactory, in view of 
the crude approximations, and confirms our interpretation of this finite size effect. 

Up to now we have assumed to be in the scaling region and therefore we have neglected 
the scaling deviations, which, at finite £o, generate corrections to the value of the string 
tension. By using the relation 



\p(B,R,L) 



1"( 



2R-1 



MP) 



(63) 



for L, R ^> £ , we extracted the value of the string tension for several values of f3. The results 
for the dimensionless quantity iV<7o£o are reported in Table | and plotted in Fig. 13 versus 
l/£o- They converge to the continuum value tt, and, as expected, the scaling violations are 
approximatively proportional to £^ 2 . From the square Creutz ratios (i.e. those with R = T) 
we obtained the same results via the relation 



xdP,R,R,L) 



MP) 



(64) 



for L, R > 



V. SIMULATIONS 

In order to check the predictive power of the 1/N expansion approach to finite-size 
scaling and to obtain results for small values of N where the expansion is not expected 
to give sensible information, we proceeded to Monte Carlo simulations of CP 7 ^ 1 models 
for a wide spectrum of values of N. In the following we will show numerical results for 
N = 2,4, 10,21,41, 100. Data for iV = 2 and iV = 4 were already presented in Refs. [|.|. 
Summaries of the results for iV = 10, 21, 41, 100 are given respectively in Tables |TJ, |HI| , [IV], 

The action @ is known to enjoy precocious scaling, and the chosen values of £ are ex- 
pected to be consistent with scaling violations 0(1%) for all the observables derived from the 
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two point function Gp In most of our simulations we employed algorithms consisting 

in efficient mixtures of over-heat bath and microcanonical algorithms. The detailed descrip- 
tion of this simulation algorithm with a discussion of its dynamical features is contained in 
Ref. JJ. Data relative to the largest lattices for N = 21 and N = 41 are taken from Ref. || 
where the simulated tempering algorithm ]TIJ was employed. 

By varying the lattice size at fixed (3 we managed to reconstruct the finite size functions 
/W and /®. They were obtained by approximating infinite lattice quantities with the 
corresponding values measured on the largest lattice available for each /3 and N (of course 
after checking the stability of the results on large lattices). The results relative to /to and 
are plotted respectively in Figs. |T^ and [15] for iV = 100, in Figs. [16] and [T7| for iV = 41, 
in Figs. [18| and [19] for N — 21, in Figs. ^ and [21] for iV = 10. In these figures the dashed 
lines show the large- A r predictions (|3T|), and the continuous lines represent the inclusion of 
the 1/N correction evaluated in Sec. |T|. 

In the scaling region the finite size scaling functions must be universal, that is inde- 
pendent of (3 and of the lattice formulation, in that they should reproduce the continuum 
physics in a periodic box. Data for iV = 10, 21, 41 coming from different f3 values show this 
universal behavior. 

A summary of all the most significant results is presented in Figs. E2] and E3] where all 
different values of N are brought together, including data for N = 2 and N = 4, in order to 
see the general trends. 

As anticipated in Sec. |T|, two physically distinct regions exist. In the small volume 
domain the 1/N expansion (as well as the perturbative expansion) is predictive for values of 
N as small as 10. The vanishing of fi in a neighbourhood of L/£ ~ 3 reflects itself in the 
approximate universality of the corresponding finite size curves in the same region for all 
N > 4. The large volume behavior shows a peculiar disagreement with the 1/N expansion: 
convergence to the infinite volume prediction is much faster than expected even for rather 
large values of N and moreover, it follows a pattern that does not seem to be amenable to 
analytic 1/N effects. 

Actually there is a plausible interpretation of this phenomenon, as we shall try to show. 
One must not forget that the physical spectrum of CP^ 1 models is characterized by effects 
of confinement, that in turn leads to the appearance of two-particle bound states, dominat- 
ing the large distance behavior of the function Gp. These bound states have binding energies 
order (6n/N) 2 ^ 3 that are nonanalytic in 1/N, and physical extension (effective radius) order 
(N/6%) 1 / 3 . Therefore they will dominate the finite size effects on large lattices (L ~ iV 1/,3 £), 
while in the intermediate region the finite size functions must interpolate between the ana- 
lytic (small volume) and the nonanalytic behavior. We may therefore expect these functions, 
for intermediate values of N, to show a (roughly) universal behavior for intermediate and 
large values of the effective variable L/ (^iV 1 / 3 ). This effect is shown in Figs. ^4| and [25], 
where this kind of universality is shown to apply to values 4 < N < 41. There is necessarily 
a saturation effect, because otherwise the convergence to infinite volume results would be 
further delayed for larger values of N, in contrast with the N —* oo prediction. Indeed at 
N = 100 we see a significant departure from the ^iV 1 / 3 scaling and an improved agreement 
with the 1/N expansion even in the large volume regime. 

At the opposite end of our analysis, when JV = 2 we find that none of the previous 
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considerations apply, but for the obvious agreement at small L/£ due to the predictivity 
domain of standard perturbation theory (that commutes with the 1/N expansion). CP 1 is 
equivalent to 0(3) non linear a model, and the finite size curves are actually showing the 
typical behavior of single-particle correlations in a finite volume (free particle in a box). 
We stress however that this different qualitative pattern seems to come out without any 
quantitative discontinuity in the dependence on N, that appears to be monotonic and smooth 
in the whole range of the variable. 

The infinite lattice correlation lengths obtained from our Monte Carlo simulations can 
be compared with those predicted by the 1/N expansion. From Eqs. (5I].3B|) we can get 
£([3) up to 0(1/ N). For N = 100 and (3 = 0.54 we find f = 1.9106, to be compared with the 
corresponding Monte Carlo result on the largest available lattice: £ = 1.922(9) (see Table 
|V]). In the continuum the 1/N expansion of £ leads to 0: 



8^ 



6.7033 ^ 



(65) 



where is the A parameter of the lattice action (||). Assuming asymptotic scaling, accord- 
ing to the two loop formula 



(66) 



we can find the correlation length as a function of (3: = (^Ai,)g(j3). For A" = 100 and 
(3 = 0.54 this would give £ = 1.955, consistent with 0(1/ N (3) deviations from asymptotic 
scaling. 

The string tension is easily extracted by measuring the square Creutz ratios Xc(R, R)\ 
indeed for L,R^> £ we should have 



9 {C) (R) 



xc(R,R) 



2R-1 



(67) 



As shown in Sec. [TIT], in practice the factor 



2R-1 



allows us to perform the measure- 



ment without requiring L ^> R. The string tension has been already measured for A^ = 4, 
N = 10 0, AT = 21, and A" = 41 ||, showing for the last two values of A" agreement 
with the large- A" prediction a£ 2 = tt/N. We have repeated the measurement for A^ = 100, 
at P — 0.54 on a L = 48 lattice, corresponding to L/£ ~ 25. In Fig. ^6] is plotted 
versus the physical distance r = R/£; the continuous and the dashed lines show the large- A" 
predictions obtained from Eq. fl43|) at the same bare coupling /3 = 0.54, respectively for an 
infinite and a L = 48 lattice. The analysis done in Sec. [IV] and the small differences between 
the infinite and the L = 48 lattice predictions show that L/£ ~ 25 should be sufficient to 
determine the string tension through the measure of g^ c \ Monte Carlo data are consistent 
with the large- A" results, which also predict the observed discrepancy of the measured string 
tension from the continuum value <r£ 2 = tt/N due to scaling violations (see Table | and Fig. 

Finally we would like to make a few remarks on the study of the topological properties 
of CP^ -1 models on the lattice. At large A^ a particularly severe form of critical slowing 
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down, exponential with respect the correlation length, has been observed when measuring 
the topological susceptibility Moreover, this phenomenon becomes worse with increasing 
N, making the simulations effectively non ergodic, already at small £. In order to overcome 
this difficulty and to perform simulations sampling correctly the topological sectors at large 
N, in Ref. || the simulated tempering method of updating |jTD[ was employed, and agreement 
with the large N result \ti 2 = 1/(2ttN) was found for N = 21 and N = 41. 

At large N we expect that the dynamically generated gauge field shows a phenomenon 



similar to that observed in the two-dimensional U(l) gauge model |TT|, that is a strong de- 



coupling in the infinite volume limit between the gaussian modes around the large- N saddle 
point, responsible for confinement, and the modes determining the topological properties. 
A manifestation of this decoupling would be the possibility to obtain good estimates of ob- 
servables not related to the topological properties even in simulations not sampling correctly 
the topological sectors ||. Relying on this decoupling property we decided to use a local 
algorithm to perform our simulations at large N. 

For N = 100 different topological sectors are separated by such a high free energy barrier 
that starting from a flat configuration, even after a long thermalization procedure, we found 
zero geometrical charge for each configuration generated in our simulations. All our N = 100 
results reported in Table [V] are effectively obtained by sampling configurations with Q g = 0. 
We also found that the geometrical charge Q g remained unchanged even when starting from 
a configuration having Q g ^ 0. As expected, the measured values of Xm, £ an d Xc(R,R) 
with R <C L become independent of Q g for large enough lattices (and for reasonable values of 
Q g ). Anyway, since the large- A saddle point configuration has trivial topology, for all finite 
lattices the correct procedure to obtain reliable results is starting from a flat configuration 
and therefore effectively sampling configurations with Q g = 0. The good agreement of the 
numerical results obtained in this way with the large- N predictions is a further support to 
the abovementioned decoupling phenomenon. 



VI. CONCLUSIONS 

The analytical and numerical results we have presented lead to some considerations whose 
relevance might go beyond the domain of CP^ -1 models. 

I) There is a perturbative domain for finite size effects, where it is possible to compute the 
finite size functions with satisfactory precision, both in the 1/N and in the standard weak 
coupling expansion. For asymptotically free theories, this domain coincides with the small 
volume regime. In the spirit of Liischer and collaborators' investigations p|,|T2||, one might 



conceive of supplementing the perturbative evaluations with numerical simulations at rather 
small volumes and matching numerical finite size effects (functions of L/£) with analytical 
results (functions of AL) in order to extract physical information (e.g. the quantity A£). 

II) In the presence of a confining potential, bound states have a physical radius that can 
be significantly bigger than the correlation length appearing in the long distance behavior. 
This new physical scale may dominate the finite size effects on volumes that are large with 
respect to the correlation length but comparable with the bound state radius itself. 

III) Phenomena related to topological properties show peculiar finite lattice behaviors. 
Actually our calculations have shown that assuming to define a topological susceptibility by 
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the standard (continuum or infinite lattice) procedure 



Xt = lim fe _o k 2 A (x) (k) , (68) 

the result would be zero on any finite lattice because the dynamically generated massless 
pole disappears on finite lattices, as expected on general theoretical grounds. As we have 
seen, this is also the reason why it is strictly speaking impossible to give a definition of the 
string tension on finite lattices. However one must not forget that the infinite volume limit 
is reached smoothly and the convergence, for sufficiently large L, becomes exponentially 
fast (water boils also in your pot, and not only in the infinite volume limit!). Therefore it 
should always be possible, with an appropriate limiting procedure, to extract infinite volume 
information from finite volume simulations. We have seen that careful handling leads to the 
possibility of "measuring" the string tension on finite lattices within the desired precision. 
Similarly, we see no obstruction to interpreting the so-called "geometrical" definitions of 
topological charge (and derived quantities) as limits of sequences of local operators @ that 
can be approximately evaluated on finite lattices, thus giving a precise meaning to apparently 
ill-defined quantities. 
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FIGURES 

FIG. 1. Large- N finite size scaling function of the magnetic susceptibility. 
FIG. 2. Large- N finite size scaling function of the correlation length. 
FIG. 3. fi^/fo^ for various values of (3. The continuus line shows the £ — > oo extrapolation. 
FIG. 4. fi^/fo^ for various values of /3. The continuus line shows the £ — > oo extrapolation. 

FIG. 5. / (P) versus i?/L at /? = 0.86. 
FIG. 6. 5 J P) versus i?/L at /? = 0.86. 
FIG. 7. 0(L/£o) at /? = 0.86. 
FIG. 8. i(/(R/L) at (3 = 0.86 for several values of L. 
FIG. 9. A^(fc,0;L) versus k£o for several values of z = L/£q. 
FIG. 10. A^(0;L) versus L/£o- The dashed line shows the function (L/£ )~ 2 - 
FIG. 11. (/>y{L/£o) for £ = 16.04. 
FIG. 12. tp Y (R/L) for £ = 16.04 and for several values of L. 

FIG. 13. The large-TV string tension on the lattice at finite £o- TV<to£o i s plotted versus l/£o- 

FIG. 14. Finite size scaling of the magnetic susceptibility for N = 100. The dashed and 
the continuus lines show respectively the leading and the next-to-leading order of the large-TV 
expansion. 

FIG. 15. Finite size scaling of the correlation length for TV = 100. The dashed and the continuus 
lines show respectively the leading and the next-to-leading order of the large- TV expansion. 
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FIG. 16. Finite size scaling of the magnetic susceptibility for N = 41. The dashed and the con- 
tinuus lines show respectively the leading and the next-to-leading order of the large- iV expansion. 



FIG. 17. Finite size scaling of the correlation length for N = 41. The dashed and the continuus 
lines show respectively the leading and the next-to-leading order of the large- N expansion. 

FIG. 18. Finite size scaling of the magnetic susceptibility for N = 21. The dashed and the con- 
tinuus lines show respectively the leading and the next-to-leading order of the large- iV expansion. 



FIG. 19. Finite size scaling of the correlation length for N = 21. The dashed and the continuus 
lines show respectively the leading and the next-to-leading order of the large- N expansion. 

FIG. 20. Finite size scaling of the magnetic susceptibility for N = 10. The dashed and the con- 
tinuus lines show respectively the leading and the next-to-leading order of the large- iV expansion. 



FIG. 21. Finite size scaling of the correlation length for N = 10. The dashed and the continuous 
lines show respectively the leading and the next-to-leading order of the large- N expansion. 

FIG. 22. Finite size scaling function of the magnetic susceptibility obtained by Monte Carlo 
simulations for several values of N. The dashed line shows the leading order of the large-iV 
expansion. 

FIG. 23. Finite size scaling function of the correlation length obtained by Monte Carlo simu- 
lations for several values of N. The dashed line shows the leading order of the large- N expansion. 



FIG. 24. Finite size scaling function of the magnetic susceptibility versus L/^TV 1 / 3 ). 
FIG. 25. Finite size scaling function of the correlation length versus L/^iV 1 / 3 ). 

FIG. 26. Rescaled Creutz ratios versus r = R/£ for N = 100. Data were taken at (3 = 0.54 and 
on a L = 48 lattice. The continuous and the dashed lines show the large- N predictions obtained at 
(3 = 0.54 respectively for an infinite and a L = 48 lattice. The dotted line is the (3 —>■ cc prediction, 
with asymptotic value tt/N. 
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TABLES 



TABLE I. Large- N string tension versus (5. 








Co 




Na il 


0.54 




2.1474 




3.2709 


0.65 




4.2859 




3.1827 


0.70 




5.8678 




3.1658 


0.75 




8.0337 




3.1557 


0.81 




11.7123 




3.1489 


0.86 




16.0350 




3.1457 


0.97 




32.0067 




3.14284 


oo 




oo 




3.14159... 


TABLE II. Summary of the simulation results for the CP 9 model. 





L 


E 


s 


Yrn 

A 111 


0.75 


18 


0.7157(3) 


3.71(2) 


19.39(11) 




21 


0.7179(2) 


3.61(2) 


18.66(10) 




24 


0.7196(2) 


3.43(2) 


17.62(8) 




27 


0.7198(2) 


3.34(2) 


17.16(8) 




30 


0.7202(1) 


3.289(15) 


16.94(5) 




42 


0.7202(1) 


3.27(3) 


16.80(5) 




48 


0.7202(1) 


3.25(3) 


16.77(4) 




54 


0.7202(1) 


3.31(3) 


16.80(4) 




60 


0.7203(1) 


3.31(4) 


16.79(4) 


0.8 


9 


0.6473(5) 


3.377(9) 


17.25(4) 




12 


0.6539(4) 


4.070(13) 


22.88(7) 




18 


0.6606(3) 


5.02(2) 


30.88(12) 




24 


0.6642(2) 


5.19(3) 


32.19(23) 




27 


0.6650(2) 


5.18(3) 


32.14(21) 




30 


0.6659(2) 


5.07(4) 


31.07(23) 




33 


0.6668(2) 


4.79(3) 


29.29(17) 




36 


0.6667(1) 


4.70(4) 


28.75(17) 




39 


0.6669(1) 


4.67(3) 


28.56(13) 




48 


0.66689(6) 


4.60(2) 


28.14(8) 




60 


0.66699(4) 


4.60(2) 


28.08(5) 


0.85 


15 


0.6145(3) 


5.221(16) 


34.20(10) 




36 


0.6213(1) 


7.21(5) 


53.63(43) 




48 


0.62196(9) 


6.63(5) 


48.73(32) 




60 


0.62216(8) 


6.48(5) 


47.41(26) 




72 


0.62214(8) 


6.44(6) 


47.18(22) 
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Summary of the simulation results for the CP 20 model. 







L 


E 




Xm 


0.65 


6 


0.7550(4) 


2.041(3) 


8.562(13) 




9 


0.7767(4) 


2.611(6) 


11.96(3) 




12 


0.7864(3) 


2.989(8) 


13.99(4) 




15 


0.7914(2) 


3.193(8) 


14.84(4) 




18 


0.7950(2) 


3.186(17) 


14.49(8) 




21 


0.7976(2) 


3.003(19) 


13.43(8) 




24 


0.7990(2) 


2.835(17) 


12.67(6) 




30 


0.7992(1) 


2.721(13) 


12.26(4) 




36 


0.7993(1) 


2.696(13) 


12.17(3) 




42 


0.7995(1) 


2.691(18) 


12.14(3) 


0.7 


6 


0.7004(6) 


2.239(4) 


9.767(17) 




9 


0.7191(4) 


2.957(7) 


14.65(3) 




12 


0.7267(3) 


3.524(7) 


18.58(4) 




30 


0.7388(2) 


4.19(6) 


21.74(25) 




36 


0.7389(2) 


3.83(3) 


19.9(1) 




48 


0.7392(2) 


3.716(25) 


19.51(8) 



TABLE IV. Summary of the simulation results for the CP model. 



p 


L 


E 




Xm 


0.57 


6 


0.8440(6) 


1.764(3) 


7.049(12) 




9 


0.8672(3) 


2.163(4) 


8.876(19) 




12 


0.8781(2) 


2.358(4) 


9.44(2) 




15 


0.8830(2) 


2.406(8) 


9.30(3) 




18 


0.8863(2) 


2.292(7) 


8.67(2) 




21 


0.8888(2) 


2.078(9) 


7.97(2) 




24 


0.8887(1) 


2.032(8) 


7.85(15) 




33 


0.8889(1) 


2.004(7) 


7.755(9) 




42 


0.8890(1) 


2.011(10) 


7.756(8) 


0.6 


6 


0.8007(6) 


1.897(4) 


7.853(16) 




9 


0.8235(5) 


2.398(6) 


10.51(3) 




12 


0.8340(3) 


2.701(6) 


11.84(3) 




15 


0.8385(2) 


2.864(7) 


12.25(3) 




18 


0.8415(2) 


2.902(10) 


12.07(4) 




42 


0.8454(1) 


2.431(11) 


10.105(15) 
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TABLE V. Summary of the simulation results for the CP model. 



P L 


E 


e 


Xm 


0.54 6 


0.8749(4) 


1.670(2) 


6.590(9) 


9 


0.8988(2) 


2.037(4) 


8.11(2) 


12 


0.9088(2) 


2.221(5) 


8.53(2) 


15 


0.9135(2) 


2.263(5) 


8.344(16) 


18 


0.9164(1 J 


z.zlU(6 J 


r .9oo(l ( ) 


21 


0.9177(1) 


2.111(5) 


7.541(12) 


24 


0.9188(2) 


2.011(7) 


7.278(13) 


36 


0.91899(7) 


1.916(8) 


7.078(7) 


42 


0.91895(5) 


1.916(8) 


7.072(5) 


48 


0.91893(7) 


1.922(9) 


7.069(6) 
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